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Abstract 

Inspired by the superblock method of White, we introduce a simple modification of 
the standard Renormalization Group (RG) technique for the study of quantum lattice 
systems. Our method which takes into account the effect of Boundary Conditions(BC), 
may be regarded as a simple way for obtaining first estimates of many properties of 
quantum lattice systems. By applying this method to the 1-dimensional free and inter- 
acting fermion system, we obtain the ground state energy with much higher accuracy 
than the standard RG. We also calculate the density-density correlation function in 
the free- fermion case which shows good agreement with the exact result. 
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1 Introduction 

In spite of the enormous recent interest in correlated fermion systems, it still remains a chal- 
lenging task to solve the interacting many fermion problem. Due to the difficulty of obtaining 
exact solutions, one often looks for a reasonably good approximation scheme for solving such 
problems. However each approximation scheme has its own limits and draw-backs. Mean- 
field methods neglect the effects of fluctuations to a great extent; the perturbative methods, 
on the other hand, do not work over the whole range of the parameter space. Among the 
other methods a real space renormalization group (RG) technique [1-7], often referred to as 
the block RG, seems to be promising in this respect. This relatively simple method includes 
some effects of fluctuations and also works over the whole regime of the parameter space. But 
the approximation involved in this scheme affects sometimes seriously the renormalization of 
the parameters (especially the off-diagonal ones like the single particle hopping amplitude). 
This, in turn, restricts the reliability of the numerical results produced by this method. 

As far as the calculation of ground state energy is concerned, the main difficulty of the 
method is the fact that by fixing a particular Boundary Conditions (BC) on a block one may 
lose a number of states which contribute to the ground state of the whole lattice. This point 
is clearly highlighted in ID Tight Binding model in which the standard RG fails for some 
type of BCs [8,9,10]. It must be noticed that this difficulty is not removed by increasing the 
size of blocks. 

In 1992 White [11] succeeded in greatly improving the efficiency of the RG method by 
inventing the Density Matrix Renormalization Group (DMRG) scheme. In this approach, 
one embeds each block into a larger(super) block and considers the block as a quantum 
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system in interaction with a reservoir (the rest of the superblock). The block can then be 
described by a reduced density matrix, whose eigenkets with the highest eigenvalues are 
used to construct the embedding operator (T) (see [12] and references therein), which is the 
mathematical artifact which reduces the dimension of the Hilbert space in each step of RG. 
Although DMRG gives very accurate results compared with previous simple RG schemes, 
its practical implementation is much more difficult and time consuming. 

Regarding the fermion systems an improved scheme to the standard RG was also proposed 
for free-fermion systems [13] which converges to the exact results of the ground state energy 
faster than the standard RG, but was not applied to interacting systems. 

In this paper we apply a modification of the standard RG scheme to a one dimensional 
interacting fermion system and obtain the ground state of the system with much higher 
accuracy compared with the standard RG [14]. For free-fermion model, our prescription 
yields much closer results to the exact values compared with those obtained from the standard 
RG and those obtained in Ref.[13], we also obtain the density- density correlation function, 
which is in a rather good agreement with the exact one. 

Our method, although does not yield as much accuracy as in DMRG, is much simpler 
practically, so that one can easily implement it on a personal computer. The required time 
is less than one second for a Pentium machine. The results are compared in table-1, table-2 
and Fig.3. 
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2 The Model 

We will consider the following ID Hamiltonian, for spinless fermions 

H = tY^{c\c i+ i + c\ +l Ci) + GY,{nin i+1 ) - ii^rii , (1) 

i i i 

where cj and c$ are the usual creation and annihilation fermion operators on the i-th site and 
the fermion number operator. When G = (free fermion case), t represents the 
hopping parameter and li is the site energy while G is a positive electron-electron repulsion 
between electrons on neighbouring sites. We will consider the electron-hole symmetric case 
where there is an average of | electron per sites (half-filled case) . In this case the Fermi energy 
Li is conveniently taken equal to G and, adding a constant, the Hamiltonian is rewritten as 
follows 

H = t£(cfe+i + c\ +lCi ) + Gj2(n t - \){n i+l - \) . (2) 

i i 

It is clear that the two terms in (2) play opposite roles. We can distinguish two different 
regimes when £ < 1 the system is conducting and < rij >= \ on each site and (ii) when 
Y t > 1 the system organizes itself in two sub-lattices A and B such that < n; >~ 0,i E A 
and < rii >~ l,i £ B. 

Our aim is to calculate the ground state energy of this fermionic system in the conducting 
regime using a modified RG prescription. In section (3) we apply the method to free-fermions 
(G = 0) and in section (4) we will consider the interacting case (G ^ 0). 
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3 Modified RG for Free-Fermions 

In the free-fermion case the Hamiltonian will reduce to the following simple form 

H = ^(4 c m + 4+iCi) , (3) 

i 

where the hopping parameter t is the only one that will be renormalized. 

The decomposition of the lattice into isolated blocks and superblocks with an odd number 
of sites per each block is shown in Fig.l. In this case the size of blocks (n B ) and superblocks 
(usb) are taken to be 3 and 5 respectively for simplicity. Therefore the block Hamiltonian 
h B and the superblock Hamiltonian h SB are 

h B = t(c\c 3 + c\c 2 + c|c 4 + c\c 3 ) , (4) 

h SB = t(c\c 2 + eld) + h B + t(c\c 5 + c\c 4 ) . (5) 

We express our prescription in the following steps : 

step-1) The superblock Hamiltonian (h SB ) is diagonalized exactly. The doubly degen- 
erate ground state is in the subspaces of particle number v = 2 and 3 , 

h SB \u = 2(3)) SB = e s B (t)\u = 2(3)) SB , (6) 

where 

= (-2.73246)t , 

is the ground state energy of the 5 sites superblock. Note that the structure of the low lying 
levels of the superblock depends strongly on the parity of nss- For usb even, the ground 
state is a singlet corresponding to the ground state of the subspace v = ^p-, while for n SB 
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odd the ground state is a doublet corresponding to the electron-hole degenerate ground state 
of the two subspaces v = n?I ^ ±1 . Here it is essential to restrict ourselves to the case usb 
odd, if we want, by keeping only two levels that the Hamiltonian conserves its initial form. 

step-2) In this step in view of the half-filled property of the system we project the de- 
generate ground state of the superblock (\v = 2(3)) SB ) onto the v = 1 and 2 subspaces of the 
block Hilbert space respectively. It is this step which effectively smoothes out the sharp effect 
of the boundary conditions by immersing the block into a superblock, i.e. the projected state is 
not restricted by any particular boundary conditions of an isolated block (open, periodic,...), 
as compared with the standard RG, but takes into account to some extent the quantum fluc- 
tuation of the rest of the lattice. Moreover by emphasizing the half-filled property, our RG 
prescription leads to a form-invariant Hamiltonian. The resulting normalized state of this 
projection is named \u = 1(2)) B . 

step-3) The embedding operator is now constructed in the following form 

r/ = (|i/=l> B (0|) + (|i/ = 2> fl (l|) , (7) 

where |0) and |1) are the renamed base kets of the effective Hamiltonian Hilbert space. 

Having the form of the embedding operator, one can calculate the projection of any 
operator onto the effective Hilbert space [12], (labelling the sites in the block 1, 2, 3 from 
left to right) we will have 

Tjc jI T I = \c/ j = 1,3, 

T}c] I T I = Xc\' j = 1,3, (8) 
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where c/ and c\ act as fermion operators in the new Hilbert space, and 

A 2 = 0.47999 , 

(a similar equation for j — 2 can be obtained, but this is not required here). 

Using Fig. 2 and the labelling mentioned above for sites in each block, one writes the 
interaction between blocks in the following form 

hf,j = t(ci lClJ + c\jc 31 ) . (9) 

Then the effective Hamiltonian between the new sites is 

h e /J = T}T}h!*T I T J , (10) 

where I and J are two neighbouring blocks, while Tj and Tj are their corresponding embed- 
ding operators. This leads to 

h e /J = t'(c\'c/ + c\'c/) , (11) 

where 

t' = \H , (12) 

and c'^j) and cj^ are fermion operators for the new sites. This equation determines the 
renormalization of the coupling constant. 

step-4) The crudest estimate for the ground state energy is just the sum of ground state 
energies of all the blocks or superblocks. Thus, 

(|) (0) = \e S B {t) , (13) 



or 

(|) (0) = \e B o(t) , (14) 

where (^) < -°- ) is the ground state energy per site of the whole system in the zeroth order 
approximation. Eq.(14) in fact defines e B as 

s B o = \s S o B ■ (15) 

However this estimate clearly neglects the contribution to the energy from interaction 
between blocks. To take this missing contribution into account, we add to (^) < -°- ) the ground 
state energy of the new lattice, and obtain 

(|) (1) = \e B {t) + \e B (0 • (16) 
Iterating this procedure we finally obtain 

rp oo -i 
JV n=0 ° 

or 

l-fssr^)- (18) 

where = t and is the renormalized coupling constant after n steps of RG. In fact, the 
above argument describes the physical interpretation of the formal limiting process defined 
in [12] in which the ground state energy is obtained as 

E = lirrin^T^ HT n . 

However in the standard RG the term e B (t) is the result of diagonalization of h B , i.e., 



T^h B T = e B , whereas in our method e B is just another expression for where e^ B 



is the eigenvalue of h SB . The results for the ground state energy per site, obtained by this 
modified RG method are collected in table- 1. The results are compared both with the exact 
results and with those obtained from the standard RG [14]. It is seen that our modified RG 
scheme yields closer results to the exact one than the standard RG method. 

We would like to comment on the non- variational nature of our results. This is caused 
by our crude estimate for the ground state energy of each block in terms of that of the 
superblock (ef = |ef B ). However one should note that this behaviour is not unusual in 
RG treatment of lattice systems. For example in standard RG at second order perturbation 
of free-fermion model [2] and isotropic Heisenberg model [16], one obtains non- variational 
results. We have also obtained the ground state energy in this case for other ratios, 
namely for the ratios | , |, ^ and y|, the results are collected in table-1. Note that our 
results in these cases are above the exact values of ground state energy. This indicates 
that the non-variational result in the | case is due to the very small block and superblock 
sizes. Therefore we expect that when the sizes of block and superblock are not very small, 
our method must indeed reproduces the variational results. Our data in table-1 show more 
accuracy compared with the standard RG results for any values of the block size. These 
data also show that the result converges to the exact one by increasing the ■^ s - ratios. 

To have a better insight about the effective Hamiltonian obtained by this modified RG, 
one can look at the nature of the constructed wavefunction. In this respect we compute 
the density-density correlation function ((njn i+r )) in terms of the distance r. In order to 
do so by our RG prescription with a block size of ub = 3, we consider a chain of length 
3 N and let N go to infinity. In this calculation the ground state of the whole system |0) 
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is replaced by T\0') where |0') is the ground state of the effective Hilbert space. Then we 
compute g{r) = (nin i+r ) for any value of r. Our results are plotted in Fig. 3 and compared 
with the exact results [17]. Our results show a very good qualitative agreement with the 
exact results while the differences in the small r is due the finite size effect of the exact 
results as mentioned in ref [17]. Although this good agreement is not specially due to our 
prescription and it can be obtained from other different schemes [2] [3] [4], this shows that 
the nature of the constructed wavefunction does not change in our scheme and the obtained 
effective Hamiltonian is a good approximate one. 



4 Interacting Fermion system 

The system under consideration is an interacting fermion system with the Hamiltonian 
defined in (2). Our RG prescription in this case is the same as in the free-fermion case. 
With n B = 3 and n S B = 5 the block and superb lock Hamiltonians are (see Fig.l) 

h B = t(c\c 3 + c\c 2 + c 3 c 4 + cjc 3 ) + G{n 2 - ^)(n 3 - -) + G(n 3 - l)(n 4 - ^) , (19) 

h SB = t(c{c 2 + eld) + G(ni - -)(n 2 - -) + h B + t(c\c 5 + c\ Ci ) + G(n 4 - -)(n 5 - i) , (20) 
Following the first three steps of our RG prescription we arrive at the effective operators 

T]c jI T I = r ] c I ' j = 1,3 , 

Tjc} I T I = r ] c\' J = 1,3, (21) 

where r] is calculated to be 

"=2^' (22) 
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Here 



a = e 4 + (1 - g 2 )e 2 + ge-2 , 
(5 = 2{e-g){2e 2 + ge-2) , 
e = - t e s B (t,G) , g=^, (23) 

and as before e^ B (t, G) is the ground state energy of h SB which is some function of t and G. 
Following Fig. 2 the interaction between blocks I and J is 

h BB j = t(cl lClJ + c\ jC3I ) + G(n 3I - i)(n u - 1) . (24) 

Then the interaction between the new sites is given by 

h e /J = TjTjhfjTjTj , (25) 

where Tj and Tj are the embedding operators of blocks I and J. In this way we arrive at 

h e /J = t'(c\ Cj + c\ Cl ) + G'(( ni - l -){nj - \)) , (26) 

where t' and G' are the renormalized coupling constants given by 

t' = rjH , 

G' = -f 2 G , (27) 

and 

~>=Jtw (28) 

After performing step-4 in the preceding section the ground state energy of interacting 
fermion system will be obtained which is given in table-2. These data clearly show that 
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the modified RG method yields much better results than the standard RG method. The 
non- variational nature of these results is due to the same reason which was described in the 
free-fermion case. Moreover, these results will be variational, if one takes larger ub and n S B- 

5 Conclusions 

In this letter we have introduced a version of RG scheme and applied it to the free and 
interacting fermion systems. In this modified scheme, we can take care of the effect of 
boundary conditions in a simple way. We have calculated the ground state energy of the free 
and interacting fermion model for different values of the block size. All of these results show 
more accuracy than the standard RG and one can see the convergence of these results to the 
exact ones. In the free case we have also calculated the density-density correlation function 
whose good agreement with the exact results verifies that the renormalized Hamiltonian is 
a good approximate Hamiltonian for the low energy spectrum of the original one. Using the 
renormalization of any operator, one can also compute the correlation function for any value 
of g(— in the interacting case. This method can also be applied to spin systems [18]. 
One may hope that for a first estimate of many properties of these systems, our method 
or its improvements (or adaptation for other lattice systems) would give quite good results. 
In those cases where application of the more sophisticated DMRG method seems to be 
hopeless, like the two dimensional Heisenberg Anti-ferromagnet, we hope our method to 
give good approximate results. Work in this direction is in progress. 
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Figure Captions 



Fig.l) Block and Super Block for the lattice chain. 

Fig. 2) The decomposition of the lattice into isolated blocks and the consideration of 
neglected bonds as the effective interactions in the new Hilbert space. 

Fig. 3) Density-density correlation function g{r) versus separation r = j — i between sites 
% and j for free-fermion chain. 
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Tables 

Table- 1. Ground state energy per site of standard RG, exact results and modified RG for an 
infinite free-fermion chain for different values of block size. 



Exact result 


Standard RG 


Modified RG 


Type of SB and B 


-0.63662(^) 


n B = 3, -0.565686 


-0.65058 


n S B = 5, n B = 3 




n B = 5, -0.5854 


-0.62256 


n,sB = 7, n B = 5 




n B = 7, -0.5966 


-0.6193 


n S B = 9, n B = 7 




n B = 9,-0.6038 


-0.6196 


n S B = 11, n B = 9 




n B = 11,-0.6088 


-0.6206 


n SB = 13,71b = 11 



Table-2. Ground state energy per site of standard RG, exact results and modified RG for an 
infinite interacting fermion chain. 



*(=S) 





0.1 


0.2 


0.3 


0.4 


0.5 


I E \ 

\ tN ' standard RG 


-0.565686 


-0.584168 


-0.603266 


-0.623008 


-0.643424 


-0.664549 


( E \ 

\fN J exact 


-0.636620 


-0.657384 


-0.679129 


-0.701826 


-0.725454 


-0.750000 


( E \ 

\ ^jy ) modified RG 


-0.650583 


-0.671674 


-0.692069 


-0.714005 


-0.736275 


-0.758952 



0.6 


0.7 


0.8 


0.9 


1 


-0.686424 


-0.709098 


-0.732630 


-0.757097 


-0.782609 


-0.775452 


-0.801804 


-0.829058 


-0.857217 


-0.886294 


-0.782419 


-0.806508 


-0.830934 


-0.855876 


-0.881525 
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